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ABSTRACT 

It is conventional to calculate the probability of microlensing for a cosmologically 
distant source based on the Press-Gunn approximation that the lensing objects 
are uniformly and randomly distributed in the intervening space with a constant 
comoving density. We here investigate more realistic cosmological microlensing 
statistics by considering the strong spatial clustering of likely lensing objects 
with each other in galaxies and their association with the clumps of dark matter 
that make up the massive halos of galaxies. Both cases in which microlensing 
objects are distributed like the observed starlight in galaxies and ones in which 
the extended massive halos themselves are also composed of compact massive 
objects capable of acting as microlenses are investigated. The distribution of 
microlensing optical depth n along randomly chosen sight lines is calculated as 
is the conditional distribution of k along sight lines near one which is strongly 
microlensed. Illustrative magnification biases are also considered. These distribu- 
tions allow us to calculate both the probability that a high redshift source will be 
microlensed in the various scenarios and the likely local k (averaged over nearby 
sight lines) at which such microlensing events will occur. Our overall result is that 
the Press-Gunn approximation is a useful order-of-magnitude approximation if 
the massive halos of galaxies are made of dark compact objects but that it fails 
badly and can be qualitatively misleading in the more likely case in which only the 
ordinary stellar populations of galaxies are the dominant source of cosmological 
microlensing events. In particular, we find that microlensing by stars is limited 
to of order 1 percent of high redshift sources at any one time. Furthermore, 
even though only a small fraction of high redshift sources are multiply-imaged 
(by galaxies), it is these sources that are most likely to be microlensed by stars. 
Consequently, microlensing by stars is usually observed at k's near 1 where the 
simple isolated point mass lens approximation is not appropriate. However, if 
CDM halos are composed of condensed objects, then more than 10 percent of 
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high redshift sources are microlensed at any given time. The vast majority of 
these sources are not multiply-imaged, and have k's smaller than 0.01. 

Subject headings: gravitational lenses: microlensing - statistical lensing - dark 
matter 



1. Introduction 

The mean optical depth to strong gravitational lensing by compact objects distributed uni- 
formly through the universe was calculated by Press & Gunn (1973) nearly three decades 
ago. The resulting probabilities are often quoted when discussing the chance that a high 
redshift source will be microlensed. In this paper we compute the distributions of microlens- 
ing optical depths obtained when compact objects are distributed as star light or in the dark 
halos of galaxies, and compare the distributions obtained to Press & Gunn's uniform result. 

Cosmological microlensing was first discussed by Chang & Refsdal (1979), Gott (1981) 
and Young (1981) who pointed out that the relative motion between stars in other galaxies 
and background sources should cause variability in the observed flux, providing an avenue 
to investigate an otherwise inaccessible mass regime. This phenomenon was subsequently 
observed in the quadruple gravitational lens Q2237-I-0305 (Irwin et al. 1989; Corrigan et 
al. 1991). The idea that most quasar variability results from microlensing due to cosmo- 
logically distributed compact objects has since been championed by Hawkins (1993). While 
microlensing is most readily observed through fiux variability, there are at least three other 
characteristic signatures. Firstly; microlensing can be detected through differential magni- 
fication of emission regions having different scales. For example, in a thorough study Dal- 
canton et al. (1994) computed the effect of microlensing by randomly distributed compact 
objects on the equivalent widths of quasar lines. Prom a lack of variance in equivalent widths 
between high and low redshift samples, they concluded that stellar mass compact objects 
cannot close the universe. Secondly; microlensing of gamma ray bursts may be identified 
through observation of a burst that repeats. This phenomenon will arise where the burst 
duration is shorter than the characteristic time delay between the microlensed images. The 
available catalogs have been searched for lensed bursts, though none have identified (Marani 
et al. 1999). Using the null result, limits have been placed on the density of randomly 
distributed compact objects in various mass ranges between 10^^^ and IO^Mq. Finally; 
mirolensing broadens the observed absolute brightness distribution of standard candles. Mi- 
crolensing of high redshift type la supernovae has been discussed by Metcalf & Silk (1999) 
and Wang (1999), and latter by Mortsell, Goodbar & Bergstrom (2001) and Minty, Heavens 
& Hawkins (2001) with emphasis on the utility of forthcoming survey samples. 
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In this paper we compute the probabihty that an image of a high-redshift source is 
subject to a microlensing optical depth k,. We also find the probability that a source is 
microlensed, as well as the conditional probability distribution of k's near the lines of sight 
to strongly microlensed sources. In Sec. 2 we consider microlensing by stars, and in Sec. 3 
microlensing by dark compact objects in galactic halos. The results from these calculations 
are discussed in Sec. 4. An important objective of the paper is to quantify the suitability of 
the assumption of randomly distributed objects. We show that the assumption is reasonable 
if CDM halos are composed of compact objects acting as microlenses, but is inapphcable for 
microlensing by stars. Unless stated otherwise, we assume a filled beam cosmology having 
Q = 0.3, A = 0.7 and h — 0.65 throughout the paper. 



2. The Probability of Microlensing by Stcirs 

2.1. stellzir microlens distributions 

In this section we describe the probability of microlensing due to the stellar populations of 
galaxies. The calculations assume a universe populated by isolated (meaning no more than 
one lensing galaxy along any line of sight) elliptical and spiral galaxies having luminosities 
that are distributed according to Schechter functions of the form 

dL \lJ ■ ^ ' 

We also assume that galaxies obey the Faber- Jackson (1976) 

L ^-^^ 



and TuUy-Fisher (1977) 
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relations, where a and v^, arc the central velocity dispersions of ellipticals and rotational 
velocities of spirals respectively. The co-moving density of galaxies is assumed constant 
with values of n*e = (1-2 ± 0.21) x 10~^h^Mpc~^ for elliptical and SO galaxies and n^s = 
(1.5±0.21) X IQ-'^h^Mpc-^ for spirals (Marzke, Geller, Huchra & Corwin 1994). The velocity 
dispersion of an elliptical galaxy is taken to be cr* = 220 km sec"^ , as is the rotational 
velocity of an L* spiral {vc^ — 220 sec~^). The constants a and 7 are given the values a — —1 
and 7 = 4 for ellipticals, and a — —0.81 and 7 = 4 for spirals. 
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The light distributions of elliptical galaxies and of the bulges of spiral galaxies are well 
described by the de Vaucouleurs profile. Assuming a constant mass-to-light ratio (T), the 
surface mass density is 



3.331 1 ' >^ 



-f— V 

Ee(r) = T SolO V ^""^ y. (4) 

The profile has a characteristic radius Rq (in pc) and density Eq (in kgm~^) given approxi- 
mately by the empirical relations (Djorgovski & Davis 1987) 

1 fj.T?^ log (ct)- 1.632 

log(hRo) = (5) 



and 



1 rx- u-2\ 1 rnoaa\ , (15.06 + 5.48 - 
log (Eo/i = log (0.366) + ^ (6) 

where 

//a„ = 2.251og(i?o) + 11.23. (7) 

For computation of the profiles of bulges in spiral galaxies, we take the 1-dimensional velocity 
dispersion to be cr = Vc/ \/2. 

The surface densities of thin disks in spiral galaxies (normalized to the Milky Way) 
can be described by a Kuzmin (1956) disk with surface mass-density of the form (Keeton & 
Kochanek 1998) 



2 , 



Erf = 0.077^ (^1 + (-)^ ) kgm-', (8) 



where is the scale length of the disk which we assume scales with v L (Bartelmann 2000) 

rrf = 350of— ) . (9) 

Wc assume that the distributions of stars in ellipticals and spiral bulges have spherical 
symmetry, and that stellar disks in spirals are circular. Wc also assume that the mass to 
light ratio is constant with radius and takes the same value for spiral bulges and elliptical 
galaxies. T is therefore related to the cosmological density of stars by 



Vt^pji — p. 



""^t^C'" r*d (10) 

P*b ~l~ P*e 



where 

da-M{a), (11) 
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Pc is the critical density for a bound universe, and M{(t) is the mass of a component with 
velocity dispersion (or circular velocity) a. Recent CMB measurements (Netterfield et al. 
2001, BOOMERANG experiment) imply a cosmic baryon density of 0.019 < h'^Q:i, < 0.026. 
Accordingly, we take 0.02 as the maximum value for VL^. 

In contrast to the observed light distribution, observations of fiat rotation curves imply 
that the total mass distribution is isothermal. We assume the total mass distribution to be 
spherically symmetric for both elliptical and spiral galaxies. The surface mass density of a 
spherical singular isothermal halo with central velocity dispersion ah = -^3/2(7 at radius ^ 
is 

where G is Newtons constant. The corresponding Einstein radius is 

where c is the speed of light, Ds and Dd are the angular diameter distances of the source 
(at red-shift Zs) and lens (at red-shift za) respectively, and Dds is the lens-source angular 
diameter distance. 

In summary, we assume that the bend angle due to gravitational lensing for both ellip- 
tical and spiral galaxies is that of the singular isothermal sphere. However, we assume the 
surface mass density of microlenses in elliptical/SO galaxies to be described by the de Vau- 
couleurs profile, while in spiral galaxies the microlenses are distributed as the sum of a de 
Vaucouleurs bulge and a thin Kuzmin disk. Thus the galaxies are composed of stars (mi- 
crolenses) and smooth dark matter that sum to an isothermal mass distribution. Below 
(Sec. 2.2) we describe the calculation for microlensing probabilities in elliptical/SO galaxies 
in some detail. Then in Sec. 2.3 we discuss the changes necessary for the calculation of 
probabilities for spirals. 



2.2. elliptical/SO galaxies 

2.2.1. the distribution of k for singly-imaged sources 

The microlensing optical depth k is defined to be the ratio of surface mass density E to the 
critical density 

^-''^^D^; ^^^^ 
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and the optical depth in stars at a radius ^ is therefore 

y;„ 3.33 (i-f 

The smooth (dark) component of dimensionless surface mass density that maintains the 
isothermal mass distribution can be written 

= - «*■ (16) 

The presence of continuous matter modifies the microlensing probability due to k^, and the 
effective optical depth to microlensing due to IS 

which is not monotonic in C,- {^K,m, ^K,out} is the set of Npairs pairs of radii having the property 
that the optical depth is larger than k, at all i^^.m < C < CK,out- Extending the calculation of 
Koopmans & Wambsganss (2001) and following Turner, Ostriker & Gott (1984), we find the 
probability that a beam from a singly-imaged (where singly-imaged refers to macrolensing 
due to the galaxy) source will pass through a microlensing optical depth (as seen by the 
observer) larger than k due to an elliptical galaxy having a central velocity dispersion between 
a and a + Aa and a redshift between and + Azd- All probabilities are calculated for 
sources at a single redshift Zg. The resulting differential cross-section is found as a function 
of central velocity dispersion a, and may be written 

^(a, K) = ^Aa^^|_,(l + z,f (max(C„ ^of - max(C, ^of) , (18) 

i=l 

where ( is the unlensed impact parameter obtained from the lens equation for a singular 
isothermal sphere 

C = e« - 6, (19) 

and is obtained from the expression for lookback time from the present (Carrol, Press 
& Turner 1992) 

^ = (1 + z)-' [{1 + z)\l + nz) - za{2 + ^)A] . (20) 

Since ^(o", z^) is monotonic in k for all za, the probability that a beam from a singly- imaged 
source will pass through a microlensing optical depth between k and k -\- Ak due to a galaxy 
between redshift za and za + Aza is found from the derivative of Eqn. 18 for each a: 

i^)-:rA:jri^^^)]- (21) 



dzddK dn \dzd 
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The magnification of a lensed image has the potential to boost the observed flux of 
a source which is intrinsically faint above the detection magnitude limit mum (Gott & 
Gunn 1974; Turner 1980), resulting in a bias for highly magnified images. The magnifi- 
cation bias for an image with co-ordinate is 

R., ^ _ N{< miim + ll0g{\l^{Q\)) .... 

Since ^{k) is multi-valued, the magnification bias with respect to an image being subject to 
microlensing optical depth between k and k -\- due to a galaxy between redshift and 
Zd + with a central velocity dispersion between a and a -\- Au is 

Note that this neglects magnification bias due to microlensing. Including magnification bias 
we find the probability that a singly-imaged source will be observed through a surface mass 
density between k and k -\- Ak due to a galaxy between redshift za and za -\- A^;^ for each cr: 

■^"""M^s.Uc.):^. (24) 



dzddK dzddK 

Using Eqn. 24 and the Poisson probability for a lensing event at optical depth Kml, 
we approximate the conditional differential probability that a singly-imaged, strongly mi- 
crolensed source will have an optical depth in stars near the line of sight of between Kml 
and Kml + Akml in a galaxy between redshift Zd and Zd + ^Zd for each cr: 

d'^PML ^ ,^ _-K..rr\ d^P^ I 



dzddKML dzddn 

d'^Pn 

~ i^ML , T'^ U=KMz, for «ML < 1- (25) 
aZdCiK 

This ignores the contribution of shear. However, at low surface mass densities microlensing in 
the presence of shear is due to an ensemble of Chang- Refsdal lenses, which have similar cross- 
sections to point mass lenses. At normalized surface mass densities near 1, the Chang- Refsdal 
diamond caustic is replaced by a continuous caustic network, which has a cross-section near 
1. We therefore feel this approximation to be sufficient for the current purpose. 

Finally, the dependences on cr and Zd are integrated out of Eqns. 24 and 25 yielding 

dP, n . r dn d^P, 



dK 



, r dn d'P^ , , 
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and 

dPML r , r , dn (fPuL 



dz, da^-^—^{a). (27) 

U-f^ML Jo Jo "0" dZdCLKML 

To compute cumulative probabilities, differential probabilities are integrated from k to oo 
because ^ ^ as /t — > oo. The probability that a singly-imaged source is subject to a 
microlensing optical depth greater than k is then 

dP 

P,{> k) = J dK'^. (28) 

Similarly, the probability that a singly-imaged microlensed source will have an optical depth 
in stars near the line of sight greater than kml is 

Pml{> Kml) = [ dn'i^L-fr^- (29) 

^ and Pk{> k), as well as j^^^ and Pml{> i^ml) due to elliptical/SO galaxies for singly- 
imaged sources are plotted in the left hand panels of Figs. 1 and 2 for Zg = 2) and = 0.001, 
0.005, and 0.020. No magnification bias (i.e. — 1) was assumed for these plots, but is 
discussed in Sec. 2.6. 



2.2.2. the distribution of n for multiply-imaged sources 

The equivalent calculation for multiply-imaged (where multiply-imaged refers to macrolens- 
ing due to the galaxy) sources was made by replacing the cross-section in Eqn. 18 with 

, , ,^.rr, ^p-- { min(CL, ^0? - min(C^, ^o)' if Qn > 

|l(a,/.) = 2^Aa^^|_,(l + .,)3 nuniCL.^of + {Qnf ^Cln<0 

' ° [ (aj'-(cj' ifcL<o 

(30) 

which is the probability that a beam from a multiply-imaged source will pass through a 
microlensing optical depth smaller than k due to a galaxy between redshift z^ and Zd + Azd 
with a central velocity dispersion between a and a + Act. ^ and Pk{> k), as well as 
and Pml{> kml) due to elliptical/SO galaxies for multiply-imaged sources are plotted in 
the central panels of Figs. 1 and 2 for Zg — S and = 0.001, 0.005, and 0.020. Here no 
magnification bias was assumed, this is discussed in Sec. 2.6. 
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2.2.3. the distribution of k for all sources 

'dnZl ' -^«(^ ^^'^ Pml{> i^ml) due to elliptical/so galaxies for all sources were found 
by summing the pairs of differential and cumulative distributions described in the previous 
two subsections for singly and multiply-imaged sources. These are plotted in the right hand 
panels of Figs. 1 and 2 for = 3 and = 0.001, 0.005, and 0.020. 



2.3. spiral galatxies 

For spiral galaxies, the optical depth in stars at a radius ^ is calculated from the sum of the 
bulge and disc surface mass densities 



■'crit 



sin[i) T^crit rd\ Td ' 



where 



e. = ^/(ecos(^)f +(-i--esin(^)) , (32) 



2 



i is the inclination of the disk (i = ^ is face on), and 9 is the angle between the image position 
and the semi-major axis (subtended at the galactic center). Sq and Rq were calculated from 
Eqns. 5 and 6, with cr = 

The bias BK{a, K,i,9) (Eqn. 23) and cross-section -^{a, K,i,9) (Eqn. 18) are computed 
as a function of i and 9. Since we have assumed that the bend angle is that of an isothermal 
sphere, the inclination and orientation averaged bias and cross-section are 

2 r- r- 

Bn{zd,cr) = —J di J d9B^{zd, a,i,9) (33) 

and 



dr , , TT^ f ^ f ^ ,„dT 
dzd 



/ di / d9—{a,K,i,9). (34) 
Jo Jo dzd 



Figs. 3 and 4 show ^ and Pk{> k), as well as and Pml{> K'Ml) for singly-imaged 

(left hand panels), multiply imaged (center panels) and all (right-hand panels) sources due 
to spiral galaxies. The source redshift was Zg = "i and Vl^ = 0.001, 0.005 and 0.020. No 
magnification bias was assumed for these plots. 



2.4. variation with source redshift 



Turner, Ostrikcr & Gott (1984) showed that the probabihty of multiple imaging by isothermal 
galaxies is a sensitive function of source redshift (cx for Zg ^ 1). We therefore expect 
that microlensing is more likely for high redshift sources. Fig. 5 shows Pk{> 
and Pml{> i^ml) for microlensing by the combined elliptical and spiral galaxy populations 
for all sources at redshifts Zg — 0.3, 1.0 and 3.0. No magnification bias was assumed, and 
= 0.005. 



2.5. veiriation with cosmology 

Turner (1990) and Fukugita, Futamase & Kasai (1990) showed that in a fiat universe, the 
presence of a non-zero cosmological constant significantly increases the frequency of multiply 
imaged gravitationally lensed quasars. We therefore expect a corresponding increase in 
the probability of microlensing. Fig. 6 shows f^^, Pk{> k) and Pml{> i^ml) for 
microlensing by the combined elliptical and spiral galaxy populations for all sources assuming 
+ A = 1 and A = 0.9, 0.7 and 0.0. No magnification bias was assumed, the source red-shift 
was Zg — 3.0, and = 0.005. 



2.6. variation with magnification bias 



To compute the magnification bias, we use the number magnitude relation described by 
Kochanek (1996) 

dN 

— = A^olO"-"*, (35) 
dm 

where A^o is a normalizing constant and is the logarithmic slope. To simplify our cal- 
culations and to keep the conclusions as general as possible, we assume that this form is 
valid for all sources at redshift Zg fainter than the survey depth. As an example, for quasars, 
and based on data from Boyle, Shanks & Peterson (1988) and Hartwick & Schade (1990), 
Kockanek (1996) finds that at Zg = 3, ~ 0.27 below ~ 19. 

In the present case of an isothermal mass distribution Eqn. 22 reduces to 



B 



(36) 



Fig. 7 shows ^^777, Pk{> and Pml{> i^ml) for microlensing by the combined elliptical 



and spiral galaxy populations for all sources at redshift z. 



=3.0, and assuming = 0.005. 
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Magnification biases were calculated assuming 3 different values of am' 0.0 (unbiased, solid 
lines), 0.2 (dashed lines) and 0.5 (dot-dashed lines). 



3. The Probability of Microlensing by CDM Halos 

In this section we assume that dark matter is in the form of compact objects and calculate 
the probability of microlensing in the dark halos around galaxies. We assume that dark 
matter halos are described by the (NFW) profile of Navarro, Prenk & White (1995, 1996, 
1997). Thus both the bend angle due to gravitational lensing, and the surface mass density 
of microlenses are described by the projection of the NFW profile. The halos are assumed 
to be isolated (meaning no more than one lensing galaxy along any line of sight), and when 
combined, to account for a cosmological density floPcrit- 

Following Bullock et al. (2001), the NFW profile has the following two parameter 
functional form for space density as a function of radius r 



which yields a microlensing optical depth in compact objects of (Bartelmann 1996) 

Kco = ^ , ,2 / ( - ) ' (38) 

where 



= < 1-^tanh-^ if x<l (39) 

if X = 1 

and Ps and are the characteristic density and radius. The concentration of the halo is 
defined as 

C.ir = — , (40) 

where r^ir is the radius containing an over-density of A^j^- In fla-t cosmologies, A^j^ is 
dependent on Q and at redshift zero is approximated by (Bryan & Norman 1997) 

187r^ + 82(0-1) -39(0-1)^ 

^vir ~ ^ ■ (41) 



The virial mass of the halo M,„v is 



vir 

47r 
3 



M^ir = —A^irPcrlir^ (42) 
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and to complete the relations between parameters, the characteristic density is 



Bullock et al. (2001) studied the density profiles of a large sample of halos in a high resolution 
N-body simulation. They assumed a cosmology having Q = 0.3 and A = 0.7, and found the 
following mean relationships: 

lo«(^V3.41og(v:^V« (44) 



where Vmax is the maximum orbital velocity (occurring at r ~ 2.16rs) and 

Taking Vmax — V^cr, we find the density profile as a function of a using Eqns. 37-45. Note 
that the assumption of an un-evolving lens population is supported in this case by Bullock 
et al. who find that is quite insensitive to redshift. The NFW halo does not contain finite 
mass. We find that cylinders of radius Rmax ~ ^Rvir around all halos contain enough mass 
to account for the density Qpcrit, and therefore assume Kco to equal zero beyond Rmax- This 
assumption only affects probabilities for k < 10~^. 

The lens equation for an NFW profile (Bartelmann 1996; Li & Ostriker 2001) is 

n^C-o^s^^ (46) 

where rj is the source position (defined in the source plane), 

^(x) = l-/(x)+ln(|) (47) 

and ^ 

(Xs = ■ (48) 

The critical image radius {^crit) for multiple imaging is the solution of 

i-^ = (49) 

which yields from Eqn. 46 the critical position (rjcrit) inside which a source is multiply-imaged, 
and the resulting critical impact parameter 

^mult — -^Vcrit- (50) 
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We have repeated the calculations of the previous section for microlensing due to com- 
pact objects in dark matter halos. Only the cross-sections are described below, as the pro- 
cedure is analogous to that for microlensing by stars, with the following differences: Firstly, 
K is more simply defined by 

and decreases monotonically with ^. Because ^(/t) is single valued, = B{^^). Secondly, 
the un-lensed impact parameter is now given by 

C = ^^viQ (52) 

together with Eqn. 46. 



3.1. the distribution of k for singly- imaged sources 

The probability that a beam from a singly-imaged source will pass through a microlensing 
optical depth larger than k due to a galaxy halo with a velocity dispersion between a and 
a + Act and a redshift between Zd and Zd + Aza is 

dr dn ctt dHot . . o ^ , .9 ^2 \ 

~dzd ^ 'd^ Ho~d^^^^^''^ ^ [maxiC^rnuit) - ^muit) ■ (53) 



3.2. the distribution of k, for multiply-imaged sources 

The mapping ([^) has 1 or 3 solutions. As ^ is decreased from large values to zero, ^(0 
crosses ( = , reaches ( = —^rnuit when = ^cHt)-! and then increases to C = 0. This results 
in three cases. The probability that a beam from a multiply-imaged source will pass through 
a microlensing optical depth smaller than k due to a galaxy halo with a velocity dispersion 
between a and a -\- Ao" and a redshift between za and Zd + Azd for each a is 

dr dn CTT dHot ( iCluit - ^HC, CmuitV) forC>0 

— {a, k) = 2—Aa——^\,=:,^(l+Zd)^ < (^^^^ + max{\C\,Uuit?) for C < and C„ > 

y) for C < and < ^crit- 
(54) 



3.3. variation with redshift and magnification bias 

dnZl ' -^«('^ Pml{> K'Ml) for all sources were found from the pairs of differential 

and cumulative distributions calculated for singly and multiply-imaged sources. Fig. 8 shows 
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these for sources at redshifts Zs — 0.3, 1.0 and 3.0. No magnification bias was assumed. 

The NFW profile is flatter than the isothermal sphere in the central regions, which 
reduces its relative cross-section, but increases the magnification of multiple images. Mag- 
nification bias may therefore be more important. For a circularly symmetric lens the mag- 
nification is 

(55) 



The magnification of an image can therefore be calculated for the NFW profile from Eqn. 46. 
Fig. 9 shows Pk{> and Pml{> kml) including magnification biases calculated 

assuming 3 different values of a^: 0.0 (unbiased, sohd fines), 0.2 (dashed lines) and 0.5 
(dot-dashed lines). Plots are shown for sources at redshift Zg — 3.0. 



4. Results and Discussion 

Probability distributions have been computed for the microlcnsing optical depth k, due both 
to microlensing by stars in galaxies, and by compact objects in the dark halos around galaxies. 
Multiply-imaged sources are generally subject to larger microlensing optical depth, though 
there is some overlap in the distributions for singly and multiply imaged sources for both 
the cases of microlensing by stars and by compact objects in CDM halos. 

First we discuss some results for the probability of microlensing by stars. Probability 
distributions were constructed for both elliptical/SO and spiral galaxy populations. For 
ellipticals we find that larger cosmological stellar densities increase the probability of each 
sometimes by a large amount. For example, the probability of each k is increased by around 
an order of magnitude between = 0.001 and = 0.020, though the increase is not 
uniform over all k, due to the factor of 1/|1 — Naturally, this results in an overall increase 
in the probability of microlensing. In addition, at higher densities microlensing is more likely 
to be observed for sources with higher Kml's near the line of sight (the distribution mode 
increases with f2*). Conversely, our model of spiral galaxies retained a constant mass-to- light 
ratio in the disc, but varied the mass-to-fight ratio in the bulge. The differential probability 
for kml < 0.1 is nearly independent of ft^, suggesting that nearly all microlensing at these 
optical depths will be due to microlenses located in the disc. However there is some variation 
with for Kml > 0.1. In particular, multiply imaged sources have distributions that depend 
on O*, suggesting that most multiple images created by spiral galaxies are located in the 
bulge. 

It is rare for a source to be microlensed by stars. Between about 1 percent (for = 
0.001) and 10 percent (for = 0.020) of high redshift sources are subject to k, > 0.01 in 



- 15 - 



stars in elliptical galaxies. As a result, only a few lOths of 1 to 1 percent of high redshift 
sources will be microlensed by stars in elliptical galaxies at any one time. About 1 percent 
of high redshift sources are subject to k, > 0.01 in stars in spiral galaxies. In our model, this 
number is nearly independent of Q^,. The resulting microlensing rate for high redshift sources 
is around 1 10th of 1 percent, comparable to that in elliptical galaxies for fl^ = 0.001, but 
lower otherwise. 

The fraction of sources that are multiply-imaged by elliptical galaxies is only around 
1 percent, and the fraction multiply-imaged by spiral galaxies is a factor of ~ 10 smaller. 
However the larger /t's near the lines of sight to multiple images means that in most cases 
where microlensing by stars is observed, the source will be multiply- imaged. For spiral 
galaxies, about 70 percent of microlensing will be observed in multiple images, while for 
ellipticals the discrepancy can be much larger, particularly if is small. Note that the 
inclusion of angular and flux ratio resolution biases will lower the multiple imaging rate. 
This will result in a microlensing rate that is higher in single than multiple images in spirals, 
and in ellipticals if is large. 

The remainder of the discussion on microlensing by stars refers to combined statistics 
for ellipticals and spirals. The probability distribution for k due to stars is sensitive to the 
source redshift. Each value of k is 10-20 times more likely along a random line of sight to 
a source at red-shift 3.0 than to a source at redshift 0.3. As a result, sources at redshift 3.0 
will be microlensed at a rate more than 100 times that of sources at redshift 0.3. On the 
other hand, we find that the most likely value of k (~ 1) is insensitive to source redshift. 

To consider the effect of different cosmologies on the rate of microlensing by stars, we 
assumed a fiat universe, and varied the cosmological constant A. Changing the cosmology 
docs not change the form of the probability distribution for k,, or the value of its mode but 
docs change the overall microlensing rate. Increasing A from 0.0 to 0.9 increases the likely- 
hood of finding a given k along a random line of sight by a factor of 5-10. In addition, the 
chance that a source at redshift 3.0 is microlensed rises from ~ 0.0015 if A = to ~ 0.025 
if A = 0.9. 

Turning our attention to microlensing in CDM halos, we find that in contrast to mi- 
crolensing by stars, all sources at redshift 0.3 are on average subject to normalized surface 
mass densities of > 0.0005, and all sources at redshift 3.0 to > 0.02. If the dark 

^crit ^crit 

matter is composed of compact objects, then this results in a microlensing rate for high 
redshift sources of more than 10 percent. The NFW profile describes halos with cores that 
are fiatter than those of isothermal profiles. The rate of multiple imaging is therefore lower, 
and this combined with the more extended distribution of mass (compared with the dc Vau- 
couleurs profile) results in microlensing by compact objects in NFW halos being dominated 
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by singly rather than multiply imaged sources. In projection, galaxy halos are likely to have 
some overlap, and a source may be viewed through more than one. However, if more than 
one halo is important, the likely values of k are low and so the approximations employed in 
the calculations of the microlensing probabilities should remain valid. 

The effect of magnification bias on the probability distribution for k was explored for 
both microlensing by stars and by objects in dark halos. We assumed a powerlaw form for 
the source luminosity function with an index a.^ (Eqn. 35). Steeper forms of the luminosity 
function boost the number of multiple images relative to single images, and for microlensing 
by stars we find the curious result that the increased magnification bias reduces the number 
of images subject to large k. This can be attributed to the fact that the highest values of 
K occur for image positions inside the Einstein radius, where the image magnification drops 
substantially. In the case of microlensing by stars, the most likely value of k near the line-of 
sight is insensitive to am- However, an of 0.5 reduces the total amount of microlensing 
by a factor of two, which rcfiects the importance of the contribution of multiple images to 
the statistics of microlensing by stars. For CDM halos, where k{^) is monotonic and the 
magnification of the central image is finite, large values of k are significantly more common 
as a result magnification bias. Because we have assumed an NFW profile, which has a low 
central density compared with an isothermal sphere, there are few multiply imaged sources. 
However these sources are significantly magnified. Luminosity functions with large values 
of am therefore significantly increase the likelihood of observing a source through a region 
of large microlensing optical depth. This can be seen from the spikes in the distributions in 
Fig. 9 (for am — 0.5). In this case (source at redshift 3), the probability of a large value 
of K can be increased more than 100 times in the presence of a magnification bias due to 
a luminosity function with a large The spikes in the distribution may be understood 
as follows. In the absence of magnification bias, radial depletion of images occurs near 
the Einstein radius where the source position lies near the optical axis, and at the critical 
radius for multiple imaging where radial arcs form. At both radii the large magnifications 
associated with the depletions result in local probabilities that are very sensitive to the source 
luminosity function. Furthermore, the Einstein radius and the critical radius for multiple 
imaging are found in the region where the profile is approximately isothermal. The values 
of K near these radii are therefore nearly constant with o", resulting in the narrow peaks 
seen in Fig. 9. The increase in the total probability for microlensing due to magnification 
bias (ttm) is only about 50 percent, demonstrating that multiple images are not a significant 
contributer to the microlensing statistics for NFW halos. 

The microlensing optical depth can be computed for objects distributed randomly in co- 
moving volume (Press & Gunn 1973; Turner, Ostriker & Gott 1984). For comparison with 
our probability distributions, we have therefore computed the values of Hunif that result 



-17- 



from the redistribution of stellar and galactic halo mass uniformly in space. For uniformly 
distributed matter the differential distributions for microlensing optical depth analogous to 
those computed in this paper are delta functions at /t„„i/, and the resulting cumulative 
distributions are step functions. These distributions are plotted (light lines) in Figs. 1-9. 
In both the cases of microlensing by stars and by dark matter, the clustering results in a 
probability distribution for k, whose mode is smaller than Hunif- This is expected since in the 
limit of high clustering, all lines of sight have zero optical depth to microlensing. However, 
the distributions have significant probability spread over several orders of magnitude in k, 
including a long tail extending to large k. As a result, microlensing by stars in galaxies is 
much more likely (> 90 percent) to be observed in images that have re's near the hne of sight 
larger than Kunif (with a most likely value around 1). The coincidence of multiple imaging 
and large values of k, (which doubles the chances of observing an image with k,), as well as 
the presence of smooth (dark) matter results in a microlensing rate for stars that is about a 
factor of 3 larger than that due to a uniform microlcns distribution. In contrast, microlensing 
in galactic dark matter halos is most likely to be observed with re's near the line of sight 
that are smaller than Kunif, although the mode lies atop a very broad peak. This results in 
a fraction of microlensed sources that (depending on magnification bias) is lower by nearly 
a factor of 2 than for uniformly distributed microlenses. 



5. Summciry 

Probability distributions were computed for the microlensing optical depth re due both to 
microlensing by stars in galaxies, and, assuming dark matter to be in the form of compact ob- 
jects, by microlenses in the dark halos around galaxies. These distributions were contrasted 
with the values of re in universes having uniformly distributed microlenses. Around 1 percent 
of high-redshift sources are microlensed by stars at any one time. In keeping with the strong 
lensing rate, incidences of microlensing by stars will be higher in elliptical/SO than spiral 
galaxies. Interestingly, multiply-imaged sources dominate the stellar microlensing statistics, 
with less than 50 percent of microlensing by stars being of singly imaged (by the galaxy) 
sources. Our model suggests that in spiral galaxies, about 30 percent of microlensing will be 
in single images, and that these will be mostly located in the disk, with microlensed multiple 
images being located in the bulge. However if CDM halos are comprised of compact objects, 
then more than ~ 1 high redshift source in 10 is microlensed at any one time. In addition, 
the vast majority of these sources are not multiply-imaged. Where compact objects are 
distributed in dark galactic halos, the simple calculation of microlensing probability from 
uniformly distributed objects provides reasonable results (to within a factor of 2). Further- 
more, the typical re near lines of sight to microlensed images in dark matter halos will usually 
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be comparable to or lower than the uniform value. In contrast, for microlensing by stars the 
typical optical depth is significantly larger than in the uniform case. Even more importantly, 
microlensing by stars will usually be observed at optical depths of order 1. 

This research was supported by NSF grant AST98-02802 to ELT. JSBW acknowledges 
the support of an Australian Postgraduate Award. 
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Fig. 1. — The probability of the microlensing optical depth due to stars in elliptical/SO 
galaxies: Differential (top) and cumulative (bottom) probabilities for k are shown. Left: 
Singly-imaged sources, Center: Multiply-imaged sources, Right: All images. The source 
redshift was Zs = 3 and no magnification bias was assumed. Functions are shown for three 
values of ^2*: 0.001 (solid lines), 0.005 (dashed lines), and 0.020 (dot-dashed lines). For 
comparison, the light lines show the corresponding probabilities where Qf.* is uniformly 
distributed. 
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Fig. 2. — The conditional probability of the microlensing optical depth along lines of sight 
near microlensed images in elliptical/SO galaxies. Differential (top) and cumulative (bottom) 
probabilities for k are shown. Left: Singly-imaged sources, Center: Multiply-imaged sources, 
Right: All images. The source redshift was -2^ = 3 and no magnification bias was assumed. 
Functions are shown for three values of fl^: 0.001 (solid lines), 0.005 (dashed lines) and 0.020 
(dotted lines). For comparison, the light lines show the corresponding probabilities where 
Oe* is uniformly distributed. 



-22- 





Fig. 3. — The probability of the microlensing optical depth due to stars in spiral galaxies: 
Differential (top) and cumulative (bottom) probabilities for k are shown. Left: Singly- imaged 
sources, Center: Multiply-imaged sources, Right: All images. The source redshift was Zs = ^ 
and no magnification bias was assumed. Functions are shown for three values of fi^,: 0.001 
(solid lines), 0.005 (dashed lines) and 0.020 (dotted lines). For comparison, the light lines 
show the corresponding probabilities where ilg* is uniformly distributed. 
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Fig. 4. — The conditional probability of the microlensing optical depth along lines of sight 
near microlensed images in spiral galaxies. Differential (top) and cumulative (bottom) prob- 
abilities for K are shown. Left: Singly-imaged sources, Center: Multiply-imaged sources, 
Right: All images. The source redshift was Zs — S and no magnification bias was assumed. 
Functions are shown for three values of Q*: 0.001 (sohd lines), 0.005 (dashed lines) and 
0.020 (dot-dashed lines). For comparison, the light lines show the corresponding probabili- 
ties where fls* is uniformly distributed. 
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Fig. 5. — The dependence of the probabihty of microlensing by stars on source redshift: 
Left: Differential (top) and cumulative (bottom) probabilities for the microlensing optical 
depth K. Right: Conditional differential (top) and cumulative (bottom) probabilities for the 
K, along lines of sight near microlensed images. Q* was 0.005, A was 0.7 (Q + A = 1), and 
no magnification bias was assumed. Functions are shown for three values of z^: 0.3 (solid 
hues), 1.0 (dashed lines) and 3.0 (dot-dashed lines). For comparison, the light lines show 
the corresponding probabilities where ft^ is uniformly distributed. 
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Fig. 6. — The dependence of the probabihty of microlensing by stars on cosmology: Left: 
Differential (top) and cumulative (bottom) probabilities for the microlensing optical depth k. 
Right: Conditional differential (top) and cumulative (bottom) probabilities for the k, along 
lines of sight near microlensed images, fl^ was 0.005, the source redshift was Zg = 3.0, and 
no magnification bias was assumed. All cosmologies had Q + A = 1. Functions are shown 
for three values of A: 0.9 (sohd lines), 0.7 (dashed hnes), 0.5 (dot-dashed hnes) and 0.0 
(dotted lines). For comparison, the light hnes show the corresponding probabilities where 
is uniformly distributed. 
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Fig. 7. — The dependence of the probabihty of microlensing by stars on the magnification 
bias. Left: Differential (top) and cumulative (bottom) probabilities for the microlensing 
optical depth k. Right: Conditional differential (top) and cumulative (bottom) probabilities 
for the K, along lines of sight near microlensed images. was 0.005, and a source redshift 
of Zs = 3.0 was assumed. Functions are shown for three values of a^'- 0.0 (unbiased, solid 
lines), 0.2 (dashed lines) and 0.5 (dot-dashed lines). For comparison, the light lines show 
the corresponding probabilities where fl^ is uniformly distributed. 
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Fig. 8. — The dependence of the probabihty of microlensing by dark compact objects on 
source redshift: Left: Differential (top) and cumulative (bottom) probabilities for the mi- 
crolensing optical depth k. Right: Conditional differential (top) and cumulative (bottom) 
probabilities for the k along lines of sight near microlensed images, fico was 0.3, and no 
magnification bias was assumed. Functions are shown for three values of Zg-. 0.3 (solid 
lines), 1.0 (dashed lines) and 3.0 (dot-dashed lines). For comparison, the light lines show 
the corresponding probabilities where VLco is uniformly distributed. 
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Fig. 9. — The dependence of the probabihty of microlensing by dark compact objects on 
magnification bias: Left: Differential (top) and cumulative (bottom) probabilities for the 
microlensing optical depth k. Right: Conditional differential (top) and cumulative (bottom) 
probabilities for the k along lines of sight near microlensed images, flco was 0.3, and a source 
redshift of Zg = 3.0 was assumed. Functions are shown for three values of a^: 0.0 (unbiased, 
solid lines), 0.3 (dashed lines) and 0.5 (dot-dashed lines). For comparison, the light lines 
show the corresponding probabilities where flco is uniformly distributed. 



